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Abstract 

By postulating the relation 9 23 — 45° + n0 13 , we seek for preferable correction terms to the 
charged lepton and neutrino mass matrices starting from the Tri-Bi-Maximal mixing and discuss 
their origins. The global fits of neutrino oscillation parameters favor n = ±l/y2; it ends up 
with the relation found by Edy, Frampton and Matsuzaki some years ago in the context of a 
T' flavor model. In contrast, the results of the v disappearance mode reported by the T2K 
and Super-Kamiokande collaborations seem to prefer 77 = 0, which turns out an almost maximal 
^23- We derive a general condition for obtaining # 2 3 — 45° + ^#13 • The condition is complicated 
by the neutrino masses and CP violating phases, so that we focus on the cases of a specific 
neutrino mass spectrum and/or characteristic CP phases. Under such simplified environments, 
we find several correction terms that take very simple forms. It is also found that the obtained 
correction terms indeed arise from the spontaneous breaking of a flavor symmetry or one-loop 
radiative corrections. 
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I. INTRODUCTION 



Now that 9 l3 of the PMNS mixing matrix has been measured very precisely by the 
reactor [l] and long-baseline Q neutrino oscillation experiments, it may be said that we 
have succeeded in acquiring a clear picture of the neutrino mixing pattern. It is the Daya- 
Bay experiment that holds the record of a precise determination of 9 l3 : the vanishing 9 l3 
is now excluded at the level of 7.7 a standard deviations with an unexpectedly large central 
value of 9f 3 B ~ 8.7° [3(. Such a large 9 13 would offer a great opportunity for us to explore 
the ordering of neutrino masses, the octant of 9 23 and leptonic CP violation. 

From a theoretical side, the discovery of a large 9 l3 might make people depressed since 
that might signal the end of a paradigm of the Tri-Bi-Maximal (TBM) mixing (4], which 
predicts 9 l3 = 0°. Nevertheless, the TBM mixing may still be available as a leading order 
one in the presence of small corrections. In fact, various ways to complement the TBM 

nn q n 

mixing by perturbing the neutrino sector [5|, |6|, the charged lepton sector [7fl or both |8J, 
have recently been proposed. In the present study, we take the same stance while paying 
special attention to correlation between 9 13 and 9 23 . 

Within a framework of the TBM mixing plus small corrections, the deviations of 9 13 
and #23 from their TBM values are given by the same parameters as we shall show in 
Eqs. (fTUl) and ( fill) : thus it is likely that they are correlated with each other. Indeed, it 
often happens in some flavor models 9|. Given this fact, it is very intriguing to recall 
the relation 9 13 = \^2\A5° — 9 23 \ found by Edy, Frampton and Matsuzaki in the context 



of a T' flavor symmetry 10] (see also Ref. ll| for earlier works). We will hereafter refer 
to this relation as the EFM relation. Interestingly, the EFM relation is now in excellent 
agreement with the global fits of neutrino oscillation parameters as demonstrated in Tab. 
HI Such a simple relation is expected to help people have an insight for model-building, 
although it should not be exact!]. 

In contrast, the r esnlts of the „, disappearance .node reported by the T2K Q and 
Super-Kamiokande 15[ collaborations seem to still favor the maximal 6*23, indicating that 



only 9 13 departs from the TBM value independently of 9 23 . If this is the case, we would 
be forced to explain the stability of 9 23 while inducing an appreciable deviation for 9 13 . 
1 At least, we do not find any underlying symmetry which guarantees the EFM relation. Thus, it should 
be affected by the running of RGEs for instance. 



Data 


flbest 
#23 


flbest 
#13 


#13 


Ref. [12] 


38.4° (38.7°) 


8.9°(9.0°) 


9.3°(8.9°) 


Ref. [13] 


51.5°(50.8°) 


9.0°(9.1°) 


9.2°(8.2°) 



TABLE I. A comparison of the EFM predictions with the global fits of neutrino oscillation 
parameters in the case of normal (inverted) mass ordering: 6^ st and #23 St are the best fit values 
from Refs. [3] and [3], while 9f^ M is a prediction of the EFM relation when employing #23 St - 

This requirement could be useful information when complementing the TBM mixing. 

In view of these thoughts, we seek for preferable correction terms to the charged lepton 
and neutrino mass matrices with a guide of 

23 ~45 O + 7A 3 , (1) 

where rj = ±1/a/2 ends up with the EFM relation, while rj = corresponds to the case 
of an almost maximal 8 23 . Note that we use ~ in Eq. (CQ), because we do not believe 
this relation holds exactly and also some approximations are used in our calculations. 
After showing our definitions of mass and mixing matrices in Sec. II, we derive a general 
condition for deriving Eq. (CQ) in Sec. III. The obtained condition is somewhat complicated 
by the neutrino masses and CP violating phases, and it does not seem to provide a simple 
correction term realizable in flavor models. Thus, in Sec. IV, we restrict ourselves to 
the cases of a specific neutrino mass spectrum and/or characteristic CP violating phases 
and show several examples yielding Eq. (Tj[|). In Sec. V, we develop two possible ways 
to realize the correction terms obtained in Sec. IV by means of a flavor symmetry or 
radiative corrections. We summarize our discussions in Sec. VI. 

II. DEFINITIONS 



We begin with the SM Lagrangian augmented by an effective Ma jorana- mass-term of 
the left-handed neutrinos: 



C = C 



SM 



h.c. 



(2) 



and consider circumstances under which 6 13 and 6 2S are enforced to be 0° and 45°, re- 
spectively, by an underlying flavor physics at leading order. Until Sec. IV, we remain 
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9 12 arbitrary in order to keep our discussions as general as possible. The leading-order 
PMNS matrix is defined as U° = {Uf) ] U® = V°P°, where U% and [7° stand for the 
leading-order diagonalizing-matrices of the charged leptons and neutrinos, respectively, 
P° = Diag(e iQ / 2 , e i/3 / 2 , 1) and 



V° 



c 12 6 12 u 

-s° 12 /V2 c° 12 /V2 -1/V2 
y s° 12 /V2 c° 12 /V2 l/y/2 ) 



(3) 



a and (3 are Majorana CP violating phases, and s° 2 (or c° 2 ) denotes sin#° 2 ( or cos#° 2 ). 
We divide mass matrices into leading and correction terms in such a way that 



M f = M? + 6M„ M„ = M° + 5M, 



where M e is a mass matrix of the charged leptons. In this definition, M° and M° are 
postulated to be diagonalized by and U®, respectively, and generate U°, whereas SMg 
and SM U give rise to slight deviations from U°. In what follows, we proceed with our 
analysis based on a premise that the correction terms are sufficiently smaller than the 
leading terms. Approximating the diagonalization of M« as 



(C/° + 5U e yM e Mj(U° + SU t ) ~ Diag(m e 2 , m% m 2 ), 

we move on to the diagonal basis of M e and redefine M u as follows: 

M V -*M V = (ttf + SUrfiM} + 6M V ){U$ + 6U t )* 

~ (^W^ )* + 5t/]M°(f/°)* + {UffM° v SUl + {UN6MMT 

= M° U +SM U , (4) 

where M° = (f/°)TM°([/°)* and we have dropped several terms in the second line. We 
stress that 5M U includes corrections stemming from both the charged lepton and neutrino 
sectors. 

One of purposes of this study is to derive SM U leading to 9 23 ~ 45° + r]9 13 ; so we go into 
more detail about 5M U as well as M°. M° is diagonalized by V°, i.e. (V°yM°(V°)* = 
Diag(m5e ia , m 2 e 1 ^ , m°) = X>°, and expressed in terms of the leading-order mixing angles 



9% and neutrino masses m2: 



(M, 



v)ll 



m\e ia {c\ 2 ) 2 



0JP( B \2 



+ m%e ltJ (s 



12) i 



(m°; 



22 



(K°)l3 



v/2 



(M, 



°)23 



-[m?e-(s; 2 ) 2 + m°e^( C ? 2 ) 2 



3J) 



(5) 



where m- are real and positive. Meanwhile, 5M U is parametrized by six complex param- 



eters, such that 



8M V 



' a d e ^ 

d b f 
{e f c) 



yet only three of a • • ■ f are independent parameters of the system. In order to see that, 
let us consider the diagonal basis of M°, namely the basis of (V°)t[M° + 5M U ](V )* = 



VI + 6V V , and define 5V V as 



' X 



5V, 



X z 
Y Z 



(6) 



where X, Y and Z are complex parameters, and the diagonal entries are omitted since 
they can be absorbed into V®. In turn, going back to the basis of M° + SM U , one finds 



a = 2s° 12 c° 12 X, 



b = -s° 12 c%X+(s° 12 Y -c° 12 Z), 
c = — Si 2 c 12 X — (s 12 Y — c 12 Z), 



V2 
1 

72 



{[( 



.0 \2 
'12) 



{s° 12 ) 2 ]X - (c° 12 Y + S ° 12 Z)} 



{[(c° 12 ) 2 ~ (s° 12 ) 2 }X + (c° l2 Y + s° l2 Z)} 



f ~ S 12 C 12^- 



(7) 



Thus, throughout this paper, we will express 5M U in terms of X, Y and Z. 
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III. RELATING 13 WITH 9 23 



By considering M^mJ 2 ]: 

MM ~ M°(M°)t + M°(5Mjt + 5M,(M°)t, 
and by regarding the second and third terms as small perturbations to M®(M®Y , are 



perturbed to be 

v eZ 
v, 3 
v r3 



c 12 



'12 



+ 



m\e~ ia X + m\e^X* 

(m 2 ) 2 - (ml) 2 
m\e ia X* + mle-^X 



(ml) 2 - (m°) 2 



'125 



m\e ifi Z* + m%Z 



(m°) 2 -(m?) 2 12 (m°) 2 -K) 



(A 2 °12' 



1 

"71 
i 

72 L 



i + 
i - 



m°e ia Y* 



6 12 



(mO) 2 - (m?) 2 
m\e ta Y* + mjjy 



12 



im!]) 2 - (m^ ~ {miy - {m% 
whereas there are no corrections to the eigenvalues upto the first order expansion. From 
these elements, the perturbed tan #12, tan #23 and sin #13 are found to be 



i?) 2 12 



K) 2 - K) 2 

m° 2 e ifS Z* +m%Z 



4) 2 



0)2 ^12 



tan 6*12 
tan #23 
sin 6>i3 = 



v e2 




Kl 




v, 3 




v T3 




v e3 \ 





t 



12 



1 + 



J) 



■Re 



m\e ia X* 



m%e~ ifi X 



Si oC 



[m 



0\2 



2Re 



12^12 

m o e ia Y * + m oy q 

0\2~ S 12 



[m 



0\2 



_ (m°) 2 - (ml) 



m° 2 e il3 Z* + m\Z 



K) 2 - K) 



(A 2 °12 



m\e ia Y* + m§y m\e^Z* + 



l^mg) 2 — \m\y { m 3 )* ~ l m 2 

respectively, and it can be seen that X is mainly responsible for the deviations of #12 while 
the deviations of #23 and 6*13 are controlled by Y and Z. 

As for 9\ 3 and #23; we are particularly interested in the relation #23 — 45° + r]9i 3 with 

±1/V2 


In view of cos #13 ~ 1, this relation may be translated into 

1 



i?) 2 " 12 



+ 



AY 12 



(9) 
(10) 

(ii) 



V 



sin 9 



23 



sin [45° + 77013] ~ -j=(l + rjsm6 13 ), 



In the first paper of Rcf . 6[ , a perturbation method is also adopted, but for M v in view of V*j 
As we shall show in Sec. IV-B, we will arrive at the same conclusion. 



y 
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leading to 



\V^\^—[l + v \V e ,\}. 
Then, by substituting V e s and derived above, one arrives at 



V 



m o e ia Y * + m^Y 



(ml) 



0\2 



[m 



0\2 ^12 



+ 



m° 2 e il3 Z* + m%Z 

Si 



[m 



0\2 



K) 



0\2 °12 



Re 



m\e ia Y* + m^Y 



m 



(ml) 2 - (m°) 2 Sl2 (ml) 2 - (m° 2 ) 2 " 12 _ 
This is the condition for obtaining 6 23 — 45° +77^13 and one of our main results in this study. 
From this condition, one may be able to extract a relational expression between Y and 
Z, which enables us to know the form of SM U . Unfortunately, without any simplification, 
it does not seem the condition provide a simple SM U realizable in flavor models. Hence, 
in the next section, we will deal with the condition under simplified environments, e.g. 
vanishing Majorana phases, quasi-degenerate neutrino masses, etc.. 



%e ifi Z* + mlZ 



(12) 



IV. DEVIATIONS FROM TBM MIXING 

From now on, we concentrate on the case of s? 2 = 1/ \/3; namely V° in Eq. (jHJ) takes the 
form of the TBM mixing pattern. Moreover, we will postulate (ml) 2 — (m\) 2 = (ml) 2 — 
(m®) 2 in Eq. ( TT2|) . As mentioned just below Eq. (jSJ), the eigenvalues are only moderately 
corrected; it may be a good approximation to identify m- with the real neutrino masses 
m i . Therefore, in the light of Amf 2 ^ ^ m 23? postulating (ml) 2 — (m®) 2 = (ml) 2 — (m 2 ) 2 
is expected to be reasonable. Even with this simplification, however, Eq. (fT2l) remains 
complicated. Thus, in the following subsections, we will restrict ourselves to the cases 
of a specific neutrino mass spectrum and/or characteristic CP violating phases. For the 
purpose of reference, we depict the running of neutrino masses as a function of the lightest 
neutrino mass in Fig. [TJ In^all of our numerical calculations, the following best fit values 
and la errors from Ref. 



Am 2 3 



12] are used: 
(2.43±o;?o) x 10 ~ 3 for Normal ordering 



(2.42io;?D x 10" 3 for Inverted ordering 
Am 2 2 = (7.54^) x 10" 5 , sin 2 9 12 = 0.307^6 • (13) 
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Normal 



Inverted 




FIG. 1. The masses of neutrinos as a function of the lightest neutrino mass for the normal (left panel) 
and inverted (right panel) ordering cases. The best fit values are input for Am^ and Am^- The black, 
red and blue curves correspond to m 1 , m 2 and m 3 , respectively. 

A. n = ±l/y/2 

• Case A-I: Y = Y*, Z = Z* and a = (3 = or tt. 

Assuming that Y and Z are real parameters and that a = (3 = or n, Eq. (Tl2!) is 
solved to be 

„ k — rjy/2 



Y (14) 

with [Y - V2Z]/i] > (< 0) for the case of normal (inverted^) mass ordering. In 
deriving Eq. (|14p . we have assumed im^ + m® = ±777,2 + m 3- A validity of this 
approximation is subject to a neutrino mass spectrum. When m 1 ~ m 2 , it should 
be applicable. Furthermore, the approximation appears to be valid in a mass region 
where m 1 is much smaller than m 2 in the case of normal ordering because 77i 3 ^> m l 2 
in that region. Nevertheless, as we shall see below, the difference between m 1 and 
m 2 gives rise to slight errors making the EFM relation somewhat hazed, k = ±1 
stems from a sign ambiguity originating in the left-handed side of Eq. (Tl2|) and 
provides us with two possible solutions. For 77 = l/y/2 (—l/y/2), Eq. ( Tl4"j) results 



Note that a minus sign appears in the right-handed side of Eq. (|12p in the case of inverted mass 
ordering. 



s 



Z=2/2Y 




FIG. 2. 9 13 as a function of 9 23 for Z = (left panel) and Z = 2\p2Y (right panel) in the case of normal 
mass ordering for Case A-I. Amf 2 , Am^ and 6 12 are restricted to be within the la bounds. 

in Z = with Y > (Y < 0) and Z = 2y/2Y with F < (Y > 0) in the case of 
normal mass ordering, yielding 



/ 



8M V 



8M V 



X 

3V2 



X 

3V2 



1 

-2 -2 
-2 -2 

1 A 

-2 -2 
-2 -2 y 



7! 



/ 



V 



-1 
1 

0-1, 

-1 1 \ 

1 -1 
1 1 j 



for 



0. 



(15) 



for Z = 2V2Y, (16) 



where X is responsible for the deviations of 9 12 and in general complex. The sign 
of Y is flipped for the inverted ordering case. Note that in the case of Z = 0, 
the dependence of Eq. ( fl2l on the Majorana phase /3 disappears; therefore, /3 can 
actually take any value in this case. 

In Fig. [21 we numerically diagonalize the neutrino mass matrix Eq. (JJJ) with Eqs. 
(I15p and (TTBT) and compute 9 is as a function of #23 for the case of normal mass 
ordering, and \Y\ are varied within ~ 0.1 eV and ~ 0.003 eV (0 ~ 0.001 eV) 
for Z = (Z = 2V2Y), while m°, m® and X are suitably tuned to be consistent with 
the la constraints of AmJ 2 , Am^ and ^ 12 . In the case of Z = the EFM 

relation is slightly hazed due to the errors of an approximation of ±m\ + mjj = 



±7«2 + rrij when m\ is small. In contrast, the Z = case does not suffer from the 
errors since it does not rely on the approximation. Almost the same figures are 
obtained for the inverted ordering case, but such a haze does not show up since m 1 
is always close to m 2 . 

Case A-II: a = (3 = in the case of quasi-degenerate neutrino mass spectrum. 
Eqs. (I14p . ( !T5|) and (Tl6|) are actually valid even for complex Y and Z in the case 
of quasi-degenerate neutrino mass spectrum when a = (3 = 0. This is because the 
imaginary parts inside | • • • | and Re[- • • ] in the left- and right-handed sides of Eq. 
ffT2j) are cancelled out, and only the real parts are constrained as 

ReZ = K -vV2 ReY (17) 

with [ReY - y/2ReZ}/r] > (< 0) for the case of normal (inverted) mass ordering. 
Then, it is observed that Eq. (fT7|) is automatically satisfied once Eq. (Till) is imposed 
for complex Y and Z. 

In Figs. |3]andHJ the scatter plots of 8 l3 vs 8 23 are shown for the normal and inverted 
ordering cases. The black, red and green dots correspond to the cases of m®,^ = 0, 
0.03 and 0.1 eV, respectively, in the case of normal (inverted) ordering. X and Y 
are treated as complex parameters with arbitrary phases, while a and (3 are fixed 
at zero. It can be seen that the EFM relation becomes clear as the scale of neutrino 
masses increases. 

Case A-III: Y = Y* and Z = Z* in the case of m 1 , m 2 -C m 3 . 

Suppose m\ and m\ are negligibly small in comparison with m®, the dependences 
of Eq. ( !T2|) on the Majorana phases a and (3 would be obscured. In this case, Eqs. 
(TL4]) , f|T5|) and ( [16]) are again valid for real Y and Z without assuming specific values 
for a and (3. 

In Fig. [5j 6*i3 and 623 are numerically computed for the normal ordering case while 
fixing m\ at zero, a and [3 run from to 2tt, and \Y\ is varied from to 0.01 eV 
(0.003 eV) for the Z = (Z = 2^2Y) case. In the case of Z = 2y/2Y , the EFM 
relation is hazed like Fig. [2} but it is more significant in the present case since 
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z=o 



Z=2/2Y 




FIG. 3. 6 13 as a function of 9 23 for Z = (left panel) and Z = 2\/2Y (right panel) in the case of normal 
mass ordering for Case A-II. Amf 2; Am§ 3 and 8 12 are restricted to be within the la bounds. The black, 
red and green dots correspond to = 0, 0.03 and 0.1 eV, respectively. 



Z=0 Z=2/2Y 




ft- [deg] 6., [deg] 



FIG. 4. Legend is the same as Fig. [3l but for the inverted ordering case. 

we are looking at a neutrino mass region where the errors of an approximation of 
±mj + 777,3 = ±777° + 777.3 i s maximized. 



B. 77 = 

• Case B-I: a = (3 in the case of m 1 ~ m 2 ^ 777°. 

In a case where 777° = m® and a — /3, the right-handed side of Eq. (112j) turns out 
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z=o 



Z=2/2Y 




FIG. 5. Legend is the same as Fig. [3J but for Case A-III. 



to be 



Re[m\e ia {Y* - V2Z*) + m° 3 (Y - V2Z)}, 



and one finds that 



Z 



Y 

71 



satisfies Eq. (112j) for rj = 0, resulting in 

U i i\ 

1 -2 -2 



A' 



3^ 



1 -2 -2 



V 



-1 11 

-1 

1 Oy 



(19) 



The results of numerical calculations for the case of normal mass ordering are dis- 
played in Fig. |6j In the left panel, is varied within 0.05 ~ 0.1 eV in which 
ml = m\ is expected to be reasonable, while m\ is fixed at zero in the right panel 
in order to see how the stability of 8 23 is affected by m\ < m\. X and Y are treated 
as complex parameters, and \Y\ and a = (3 run from to 0.007 eV and to 2ir, 
respectively, in both cases. It can be seen that 8 23 scarcely departs from 45° in 
the left panel, whereas 6 23 can be 45° ± 1° in the right panel. As for the inverted 
ordering case, the assumption m® = m\ is always valid, yielding the almost the 
same figure as that in the left panel of Fig. El Hence, we refrain from showing it. 
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6„ [deg] 6., [deg] 



FIG. 6. 9 13 as a function of 23 in the cases of normal mass ordering for Case B-I. Amf 2 , Ato^ and 6 12 
are restricted to be within the la bounds. In the left panel, is varied within 0.05 ~ 0.1 eV, while it 
is fixed at zero in the right panel. 

We add a remark that Eq. ()18p is equivalent to the condition found in Ref. It 
can be seen from Eq. © that Z = Y/y/2 (with s° 13 = l/y/3) results in b = c; this 
is what observed in Ref. jg]. 



• Case B-II: -Y = Y"*, -Z = Z*, a = or tt and /3 = or tt. 

Suppose that Y and Z are pure imaginary and that a and f3 are equal to or tt, 
the right-handed side of Eq. fTT2|) turns out pure imaginary, and thus Eq. (JT2j) is 
satisfied for = 0. 

In Fig. [71 we plot as a function of # 23 (left panel) in the case of normal mass 
ordering for \Y\ = ~ 0.013 eV, \Z\ = ~ 0.001 eV and m? = ~ 0.1 eV. This 
case prefers a somewhat small X. In the right panel, we show the dependence of 
# 23 on \X\; one can see that the deviations of 8 23 from 45° become large as \X\ 
increases. The red line in the left panel corresponds to the case of X = 0. Figures 
for the inverted ordering case are almost the same as Fig. UJ so we refrain from 
showing them. 
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FIG. 7. 6> 13 as a function of 6* 23 (left panel) and 9 23 as a function of \X\ (right panel) in the case of 
normal mass ordering for Case B-II. Anij 2 , Am^ and 9 12 are restricted to be within the la bounds. In 
the left panel, \X\ is varied within ~ 0.001 eV, and the red line corresponds to the case of X = 0. 

We mention about the ambiguity of a form of 5M V : one can subtract terms proportional 
to M° and/or the unit matrix from SM U without any modification to mixing. For example, 

/ 



X 

3V2 



4 1 1 | 

1 -2 -2 
1 -2 -2 



can be divided into 



X 

3V2 



3 2 2 
2 -3 -1 
2 -1 -3 



» ^2 ^ 



+ 



2 
2 



d 1 A 

1 1 1 
1 1 1 



/J 



The second term does not affect the mixing, and the third term can be embedded into 

/ Q O O \ 



M°. In this sense 



X 

3V2 



d 1 1 \ 

1 -2 -2 
1 -2 -2 



and 



X 

3V2 



3 2 2 
2 -3 -1 
2 -1 -3 



are equivalent and result in the same mixing. The difference between them is the 
parametrization of SM U , and it is visible in the basis of 5T> U in Eq. (E]). On one hand, the 
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original form is translated as 



/ 



(V 



X 

3V2 



4 1 1 \ 

1 -2 -2 
1 -2 -2 



On* 



X 



( 1 o\ 

1 




as we defined. On the other hand, the subtracted version turns out 



(V 



X 

3^2 



(32 


2 \ 


2 -3 


-1 


^2 -1 


- 3 J 



0\* 



X 



V 



1 i73 
-fy 



and non-zero values are entered in the diagonal elements. Namely, only the eigenvalues 
are altered by this change. 



V. POSSIBLE REALIZATIONS OF 5M U 



We here describe two types of procedures for realizing Eqs. (fl~5l) . (fT6|) and ( |T9|) . The 
first one utilizes a flavor symmetry, while the second one is based on radiative corrections. 



As for Case B-II, it may be interesting to invoke the timeon mechanism 16[, which is a 



mechanism of spontaneous CP violation triggered by the VEV of a pseudo-scalar. 



A. Flavor symmetry approach 

We adopt the 5*4 symmetry as our flavor symmetry following the definitions given in 
Ref. [ill and suppose that the three generations of the left- and right-handed leptons 
belong to 3i of S4. Besides, in order for the S4 flavor symmetry to be spontaneously 
broken, the following scalar fields and vacuum expectation values (VEVs) are introduced: 

(<t>3 1 ) = (v:,v:,v a n + (t»x 2 x 3 ), 

(K) = v c u (20) 
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for the neutrino sector, and 

(4) = (Kf:0,o) + (tf,44), 
(0U = (O,o) + (<</ 2 y 3 ), 
(4) -(4,4), 

= (21) 

for the charged lepton sector. The subscripts beside <p u,e indicate representations under 
the S4 symmetry, but we do not specify their representations under the gauge symmetries. 
Here, we are only interested in how the £4 symmetry is broken by each representation 
and what sort of breaking terms are arisen, so that we mention nothing about the gauge 
structures and consider the effective Yukawa interactions 



A v L v L<t>A 



A 



J2 Y ^ e L e R^ + hx - 

A 



(22) 



where v L is the left-handed neutrino, e L and e R are the left- and right-handed charged 



leptons, respectively, and stands for a Yukawa coupling. For example, one can imagine 
Y lMii e R if 0ii is SU{2) L doublet. Meanwhile, (f/A 3 )LLHH{ 



^31^31)2 could be invented 



if 03 is a gauge singlet scalar, where L is the left-handed doublet lepton, H is the 



SM Higgs, / is a dimension-less coupling and A denotes a cut-off scale. In this case, 



/A corresponds to $ v 2 and f(H) 2 /A 2 = Y 2 U in Eq. f l22|) . 



Note that, in the neutrino sector, we have omitted 05l 2 and 0g 2 since I2 and 3 2 are not 
produced from ^(3i) (8) v* L (Zi). Similarly, <f)[ 2 have been omitted in the charged lepton 
sector. 

We presume that V£ b c and V* b c are sufficiently larger than the other VEVs and that 
they constitute the leading-order mass matrices: 

/ 

M° = 

v a 

-V" + v h u ■ 



2V U 

v a 



v v 

' c 



-V u 

a 



V b 
V b 



-V v 

y a 



\ 



and 



/ 



\ 






V? + Vc 

V? 2V? + W 



(23) 








.yi _ yi + yi 

a b ' c 



\ 



-V + K + V 



(24) 
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where the Yukawa couplings and are embedded into the VEVs. The same shall 
apply in the following discussions. On one hand, Mf is diagonalized by rotations of 
the right-handed fields, resulting in mj = 2V£ + V*, m° = — V£ + V£ + V* and mj = 
—V* — Vi? + V* as well as = 1. On the other hand, by comparing M° with Eq. 
(J3J), one immediately observes that M° takes the same form as M° for s® 2 = 1/V3) i- e - 
= ±[m?e iQ +m§], = |[-m?e ia +2m!^+m!]] and ^ = \[m\e ia +m° 2 e^ -m§]. Thus, 
these mass matrices lead to the TBM mixing, and we identify the above M° with M°. 
As it is well known 18[ , the mass matrices possess residual Z 2 x Z' 2 and Z 3 symmetries in 
the neutrino and charged lepton sectors, respectively. They are the symmetries ensuring 
the TBM mixing. 

The other VEVs break the residual symmetries and yield the correction term 5M V . In 
the charged lepton sector, we have 



6M, 



( d{ 4 ^ 

d{ d l 2 

\4 o d[j 



I 



+ 



2t[ -t 



t s 2t 2 



' i 



/ 



+ 



^ t 2 t{ 2t 3 J 



— Un 



\ 



Hi 



-u. 



\ 



J 



(25) 



Let us convert them from the charged lepton sector into the neutrino sector along the 
procedure outlined in Sec. II. In the limit of m° e = and |5M^| 2 = 0, M(M\ is given by 



MoM, 



( 





|m°| 2 

I [X I 










\m 



012 



/ 



(ml)*(4-ti-ui) 



2Re[m°i 



u 



[m 



0\* 



(d{ + 2t{ 



+ 



m° T (d\ 






- 4 + u{) 








(4 



[m 



(d 



m° T (di + 2ti)* -2Re[m° T (t{] 



ml{d\ + 


4 + 
+ 24) 

2Re\m°< 



2t$ 



\ 



J 
\ 



By regarding the second and third terms as small perturbations to the first term, 5U t is 
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found to be 



8U, 



-t%—u 



"1 



4+uj \ 



m° 
dl+2t 



nP 



-tS+uf 



ii+24 







/ 



where we have dropped terms proportional to m^/{m T ) in the 23 and 32 elements. Now 
we are ready to calculate 5U\M u {U t )* + {U^M^SU^, but the expression is long and very 
complicated. Hence, we here suppose CP invariance in the charged lepton sector; then it 
is simplified to be 



5UlM°MY + {U^Ml5Ut 



I 



12 



my 



-2(-K l 

-V" + V?) - 



vn 

-V v 



\ 



+ 



K 23 

m2 



~(-Va+ V c) 


■Va + V?) -2( 



* 

-V" 




2(-V? + Vi 



:-v: + vt 



* 




\ 



+- 



V 
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my 



V 



-2{-v:+v b n 

-v n v + vn-{-vr 



2(-K l 

* 




'Va + V? 



-V v 

y a 



\ 



J 



(26) 



where V 



12 



4-4-4, v. 



2) "23 



4 + 24 and Vf 3 = d{ 



Un 



{t[ and u[ only contribute 



to the charged lepton masses.) 

Corrections from the neutrino sector are simply given by 

/ 

u u x u 2 



2t\ t\ t y 2 



tg 2^2 



' 1 



(27) 



< ^ 
\d% d[J y 

In summary, <5M^ consists of Eqs. (1261) and (12"T|) in this scenario. For instance, in a 
case where a — (3, the first term of Eq. (1261) turns out to be the second term of Eq. (1191) 
since — + is much smaller than —V// + V/f in this case. Similarly, the second term 
of Eq. (fl6|) can be obtained by combining the second and third terms of Eq. (|26|) when 
Vf 3 = 2V| 3 . As for Eq. (ETJ), (^,t£,^) oc (2, -1, -1) and «,^) oc (1, -1) yield the first 
and second terms of Eq. (IT5|) . respectively. 
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B. Radiative corrections 



We consider the combination of tree-level and one-loop neutrino masses. The tree- 
level neutrino mass matrix is assumed to take the form of M° in the diagonal basis of the 
charged lepton mass matrix, leading to the TBM mixing at tree level. For this setup, we 
insert two types of one-loop radiative corrections and show that they may be useful to 
realize a preferable 8M V . 

The first example is proposed in Ref. 



19], in which we have 



(28) 



at one-loop level, where D e = Diag(m e , m^, m T ), v ew = 174 GeV and J 1o °p is a dimension- 
less parameter including some loop factors. The details of / 1o °p depend on model details, 
so we remain it arbitrary in this study. If a = and (3 = it in the case of quasi-degenerate 
mass spectrum, Eq. (128]) turns out to be 



3 



m, 



( 



\ 



0-2 
-2 2 -2 
0-2 



\ 



+ ml 



( 



\ 



0-2 
0-2 
-2 -2 2 



\ 



(29) 



where m° u 
by adding 



mi 



777.3, an( i we have omitted the term proportional to m\. Then, 



— (m M + m r ) 



/ 



V 



-1 -1 \ 
-1 1 -2 
-1 -2 1 



which can be absorbed into M°, one obtains the second term of Eq. (115]) : 



/ 



V 



-1 1 \ 
-1 1 

1 -1 j 



(30) 



This example may be applicable to Case A-Ic 



Recall that there is no constraint on j3 in the case of Z = 0. 
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The second example is the so-called Zee-model given in Ref. 
/ 



5M„ 



fe^l-^l) /erK-We) 



201 ] . which induces 
\ 



f er (m 2 T - m 2 e ) f^{m 2 T -ml) 

a 
a c 
b c J 







1 

r 



/ 



at one-loop level. f a/3 are Yukawa couplings, and T contains loop factors and possesses 
mass-dimension one. Similar to the first example, by adding and subtracting terms pro- 
portional to 

/ n n n \ 



/ 



1 1| 

1 1 
1 1 y 




1 -1 
-1 1 J 



( 



1 
1 
1 



\ 



one arrives at Eq. (JT9J): 





U i 


i \ 




{ ° 


1 






8M V = X 


1 -2 


-2 


+ Y 


1 








(31) 




I 1 - 2 


- 2 ) 




I- 1 





) 





with X = (a + b - 2c)/18 and Y — (a — 6)/2. 



VI. SUMMARY 



The TBM mixing may still be applicable as a leading order one in the presence of 
small corrections, although the exact TBM mixing was ruled out by the discovery of a 
non-zero 6 13 . If this direction is adopted, one immediately faces with a difficult challenge 
of revealing structures and origins of the corrections. In this work, employing Eq. ([T]) as 
a guide, we seek for preferable correction terms to the charged lepton and neutrino mass 
matrices and discuss their origins. The latest global fits of neutrino oscillation parameters 
point to T) = ±l/v2, whereas the results of the v disappearance mode seem to prefer 
t] = 0. We have succeeded in deriving the general condition Eq. ( fl2l) for obtaining Eq. (00), 
but it has also been found that the condition is complicated by the neutrino masses and 

20 



CP violating phases. Hence, in order to find out simple and realizable correction terms, 
we concentrate on a specific neutrino mass spectrum, e.g. the quasi-degenerate spectrum, 
and/or characteristic CP violating phases, e.g. vanishing Majorana phases. Under such 
simplified environments, we arrive at Eqs. ( 1L5|) and ffl~6]) for the case of 77 = ±l/\/2, 
while Eq. ( TT9|) for 77 = 0. Furthermore, we develop several ways to realize the obtained 
correction terms by utilizing the spontaneous breaking of an S4 flavor symmetry or one- 
loop radiative corrections. 
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